Exact Bogoliubov-de Gennes solutions for grey soliton backgrounds 
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We derive and discuss the complete set of exact solutions to the one-dimensional Bogoliubov-de 
Gennes equations for small amplitude excitations around general grey soliton solutions to the cubic 
nonlinear Schrodinger equation. Our results extend the previously known case of the motionless dark 
soliton background. We derive our non-zero frequency solutions using a variant of the factorization 
method for Schrodinger equations with refiectionless potentials. We also discuss the zero mode 
solutions at length. 
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I. INTRODUCTION 

Solitons are a special class of solutions to a special 
class of nonlinear wave equations and appear as such in 
many different physical contexts. They are remarkable 
because they are, in a qualitative way, typical of a phys- 
ical phenomenon so fundamental and ubiquitous that it 
is often entirely overlooked: The few parameters of a 
soliton, such as its position and velocity, stand out, as 
natural collective coordinates, from the infinitely many 
other degrees of freedom of the field in which the soliton 
propagates. Whenever physicists model something real 
as a finite idealized 'system', we rely implicitly on the fact 
that a few particular degrees of freedom often do stand 
out in some way from among many. And whenever we 
speak of an open system, subject to noise or dissipation 
or even quantum decoherence, we are stating implicitly 
that the environment is somehow less prominent than the 
system. Prominence or its absence are thus very often 
taken for granted, which is more easily accepted than ex- 
plained. It would be convenient to leave the issue to the 
philosophers, if it did not have such important physical 
consequences. It is worth looking at even simple mod- 
els for prominence within physics. Solitons are a special 
case in which this general but elusive phenomenon can 
be clearly examined. 

The first step towards understanding the dynamics of 
solitons, as preferred collective coordinates among in- 
finitely many other field modes, is of course to identify 
the soliton solutions themselves. This has been done in 
very many cases, and the results constitute much of the 
extensive literature on solitons [lU3|. The next step, how- 
ever, is to learn about the solutions to the field equations 
within a configuration space neighborhood of the soli- 
ton, by solving the equations linearized about the soli- 
ton background. Far fewer explicit analytical results are 
known for this higher level question. The contribution 
of this paper is to provide them for a nontrivial one- 
parameter family of solitons studied in nonlinear optics. 



in hydrodynamics, and in Bose-condensed dilute gases - 
the grey solitons of the cubic nonlinear Schrodinger equa- 
tion (NLSE) in 1+1 dimensions (one dimension of space, 
one of time) [Sj]. 

Since grey solitons have been observed in quasi-one- 
dimensional Bose-Einstein condensates, where their mo- 
tion has been confirmed [H-llI to follow quite well the pre- 
dictions of the NLSE in its role as the Gross-Pitaevskii 
mean field theory [9] , there have been several recent theo- 
retical studies which investigate solitons as more than iso- 
lated NLSE solutions. Various numerical methods have 
been applied to evolve many-body quantum states of 
a dilute Bose gas, and identified discrepancies between 
the quantum evolution and its mean field approximation 
given by the NLSE. This concerns the dark soliton dy- 
namics in a lattice [lol-[l3. or in a combined harmonic 
and lattice potential IsT 14[ also at finite temperatures. 
In (T5I4T7I the dissipation of the soliton due to thermal 
fluctuations explained the experimentally observed in- 
crease of the soliton oscillation in the harmonic trap [l^ ■ 
Lately, [lOl] used an effective Lagrangian for the dark soli- 
ton to estimate its decay and flnite lifetime due to scat- 
tering on background fluctuations, which are treated in 
a Luttinger liquid approach. 

The complementary work reported in this paper is in 
one respect much less ambitious than these quantum 
studies, and in another much more. We work entirely 
within a linearized theory in which there is not even any 
difference between quantum and classical dynamics, and 
so we ignore the entire issue of quantum corrections to 
the NLSE. On the other hand, though, our results are an- 
alytical, explicit, exact, and Hamiltonian. We therefore 
provide tools which may in future be used to understand 
the quantum dynamics that can be revealed, or at least 
suggested, by numerical many-body calculations. 

The linearization of the NLSE about a stationary 
background is known in the dilute gas context as the 
Bogoliubov-de Gennes equation (BdGE). BdGE solu- 
tions are already known for the so-called 'bright soliton' 
backgrounds [20'| that are found for the case called 'self- 



focusing' in optics, and 'attractive' in gases [21[. Bright 
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solitons are in a sense a trivial limit of solitons as col- 
lective degrees of freedom, however, because away from 
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the soliton, the field amphtude falls off rapidly to zero. 
The bright soliton's position is thus simply a sort of cen- 
ter of mass for the entire field, and its motion is simply 
a Galilean boost of the entire system. A direct conse- 
quence is the gap in the Bogoliubov excitation spectrum, 
which simplifies the analysis of the soliton lifetime, for 
example, in contrast to the ungapped spectrum of the 
repulsively interacting gas p2j . 
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FIG. 1: A grey soliton located at a; = 0. Shown is l*!/^ for 
periodic boundary conditions. The precise depth and breadth 
of the soliton 'notch' are determined by the speed with which 
it moves. The Bogoliubov-de Gennes problem concerns the 
infinitely many modes of perturbations in 'i!{x,t) around this 
background. 

With oppositely signed nonlinearity coefficent, the cu- 
bic NLSE is instead 'self-defocusing' in optics, and repre- 
sents repulsively interacting Bose-Einstein condensates. 
In this case it has 'grey soliton' solutions, in which 
the field amplitude is locally reduced below its non-zero 
asymptotic value; there is a 'notch' in the wave function 
modulus l^*!, as shown in Fig. [U with a concomitant lo- 
calized twist in the phase of ^f. The grey soliton position 
is thus distinct from the background center of mass, and 
its velocity relative to the field at infinity is an indepen- 
dent degree of freedom from the global velocity of the 
entire system. BdGE solutions have recently been pub- 
lished for the strictly 'dark' soliton, which does not move 
relative to infinity p3l [2^ . In this paper we present the 
BdGE solutions for grey solitons of all velocities in detail 
and at length, for periodic boundary conditions as well 
as on the infinite line. We exhibit a continuous spec- 
trum of BdGE normal modes, plus discrete zero modes 
(normal modes of zero natural frequency) that are asso- 
ciated with perturbations of the soliton motion itself. We 
demonstrate the completeness of our modes as a basis for 
all BdGE solutions. 

With the Bose-Einstein condensate foremost in mind 
as the physical realization of our equations, we analyze 
the BdGE as the canonical equations of motion of a 
Hamiltonian system, and discuss the amplitudes of our 
normal modes as canonical coordinates. This terminol- 
ogy should properly be modified for optical or hydro- 
dynamic realizations of the NLSE and its linearization. 



since in these cases the variable we denote as t is actu- 
ally spatial distance, and so the functional we call the 
Hamiltonian generates translation instead of time evolu- 
tion. The periodic case that we refer to as a 'ring ge- 
ometry' would also rather apply, in optical contexts, to 
solutions in an unbounded spatial line but with a given 
temporal period. All our results remain formally valid, 
however, for any version of the BdGE equations. It is 
likewise unimportant for the present paper whether the 
dynamical system represented by the BdGE is quantum 
or classical. We are simply separating linear equations 
of motion by finding normal modes, and so our results 
are the same whether the amplitudes of those modes are 
operators or c-numbers. 

The paper is organized as follows: Throughout the pa- 
per we will present results for both infinite and periodic 
systems, considering the two cases always in parallel. In 
order to provide a self-contained and pedagogical discus- 
sion of the problem, the first half of the paper is a review 
of standard material. In Section |TT] we review the NLSE 
as a Hamiltonian field theory, and in Section IIIII we de- 
scribe the grey soliton solutions themselves. In Section 
IIVI we review linearization around a background NLSE 
solution and present the BdGE as a linear Hamiltonian 
system. In Sections IVland rvTl we finally present our main 
results, the explicit BdG eigenmodes for grey solitons of 
arbitrary velocity. A reader already familiar with both 
grey solitons and BdG theory might begin reading here. 
Section IVII is devoted to the zero modes, which are of 
special interest and subtlety. Even readers quite familiar 
with BdG problems may wish to consult Section IIVI be- 
fore Section lVll in order to review the less commonly used 
canonical representation we adopt in this paper. Sec- 
tion |VlIl discusses our results and suggests further lines 
of study. A series of appendices treats technical points 
skipped in the main text. 

II. HAMILTONIAN 

We begin with the Hamiltonian functional of a complex 
field ^(a;,t) in one spatial dimension: 

Ho^Jdx . (1) 

As dynamical variables in Hamiltonian field theory, ^ 
and 5** are symplectic co-ordinates (25j : their real and 
imaginary parts correspond respectively to the usual 
canonical position and momentum fields. Consequently 
the canonical field equation for ^ derived from ([1]) is the 
nonlinear Schrodinger equation 

This equation may also be obtained as the Heisenberg 
equation of motion for the second-quantized destruction 
operator for a dilute Bose gas with repulsive particle in- 
teraction in the s-wave scattering limit, in appropriately 
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re-scaled dimensionless variables, with the quantum field 
^ then replaced with the c-number field 5* in Gross- 
Pitaevskii mean field theory. Nevertheless, it is worth 
appreciating that Eqn. ^ is also a perfectly good clas- 
sical Hamiltonian field equation, and appears as such in 
several physical contexts. 

Our goal is to describe the linearized dynamics gener- 
ated by Hq around a time-independent solution to the 
field equation - namely, a soliton solution. But in fact 
we can take advantage of two symmetries of Hq, and 
obtain linearized dynamics around solutions with gener- 
alized time-independence, as follows. If we define a new 
field tlj{x,t) such that 

^{x, t) = e-'(^+'"^-'''/2)V(2: - (/3 - v)t, t) (3) 

for constants /?, v and /x, we can describe any field 4* in 
terms of ip instead. 

The field equation for ip is derived straightforwardly 
from and after changing to a moving co-ordinate 
frame (a;', t) = {x — [fi — v)t, t) reads 

idMx\t) = -^dl,iP + - v)d^,^ -f IVP^ - , (4) 
where we have defined 

li + vP ~ — (5) 

to shorten the equations. Throughout this paper we will 
work in this moving frame and drop the ' to simplify 
notation. The corresponding Hamiltonian functional of 
-0, whose symplectic field equation is (jH) just as that of 
Ha is ©, is 

H = jdx Qia.vf + 

o _ X (6) 

The differences between H and Hq are the constants of 
motion associated with the symmetries of phase and po- 
sition translation, and the change from and Hq to ip 
and iJ is a canonical transformation. In terms of a Bose 
gas, the operators at issue here are the total particle num- 
ber and total linear momentum, which commute with the 
Hamiltonian. Evolution within the subspace of definite 
number and momentum may be generated equivalently 
by the quantized versions of H and Hq^ since the dis- 
crepancy between them is only an unobservable phase 
prefactor mulitplying all amplitudes in the subspace. 

The main result of this paper will be to explicitly diag- 
onalize the linearization of H around a non-trivial class 
of time- independent solutions to (|3]): grey solitons. 

III. GREY SOLITONS 

The grey soliton [ll] with position xq (constant in 
our CO- moving frame), and constant velocity /3 relative to 



the background phase gradient, is 

^q{x) = e*''e-™(^^="°Hi^ + KtanhK(x - a;o)} 

The constant overall phase 9 can be chosen to be zero 
without loss of generality. For /tjcc — xq] ^ 1, ip has the 
constant phase gradient v corresponding in the BEG con- 
text to a constant gas velocity v. We will sometimes refer 
to this solution parameter as the (constant) gas velocity 
in the following. The soliton itself is moving relative to 
our original frame slower by (3 (or faster by — /3) than the 
background gas motion. It is thus important to recognize 
that /3 and v are independent parameters. A grey soliton 
moves with the phase gradient of its background field, 
but it is also able to move independently of its surround- 
ing field, by deforming its shape. Note that the relative 
soliton speed /3 is restricted to lie within the range 

I/3|<Vm = c, (8) 

i.e. the soliton speed can never exceed the speed of sound 
c of the background gas. For the extremal values of /3, 
the grey soliton reduces to a translationally invariant tp 
with constant modulus and phase gradient. 

The squared modulus ['(/'(a;)!^ (corresponding to con- 
densate linear density, or light intensity in a nonlinear 
fiber) has a minimum at x = xq, where it drops to 

l^olL„ = (9) 

The case /3 — Q is therefore called the 'dark soliton' (since 
in fiber optics it represents a perfectly dark spot) al- 
though grey solitons with nonzero /3 are also sometimes 
referred to as dark. \tpo\'^ rises with |a; — xo\ toward its 
asymptotic value of -I- /?^. Except for the extremal 
values of /?, for which k becomes very small, grey soli- 
tons are well localized objects. The asymptotic value of 
I ■(/'op is nearly reached already by around k\x — xo\ > 2, 
and the remaining deficit -I- — [■(/'op then decays 
exponentially on the soliton length scale 1/k, see Fig. [T] 
The soliton solution can also be characterized by its 
phase behavior. Across the interval \x — xa\ < d, ip has 
a total phase difference 

Apd = ±7r - 2 arctan [ ) + 2vd (10) 

\ K tanh nd J 

for /3 ^ 0. Apart from phase change due to the gas 
velocity u, the largest jump in phase occurs over Kd ~ 2. 
The grey soliton is localized in its phase behavior as well 
as in its modulus. 

A grey soliton with larger (3 is moving faster relative 
to its background, but it is a longer and shallower dis- 
turbance of the background. For this reason it has lower 
energy H than a dark soliton with /? = 0, which does not 
move relative to its background at all. Gonsequcntly one 
may state that dark solitons are energetically unstable to 
acceleration |j9|]! 
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Periodic boundary conditions may also be imposed on 
our system; in the BEC context, for example, the gas may 
be confined in a tight toroidal trap [27l - l29j . The soliton 
solution ([7]) is in this case no longer exact, because it 
cannot be made smoothly periodic; the exact periodic 
solutions may be expressed in terms of Jacobi elliptic 
functions , and approximate solutions can be obtained 
[soj . In this paper we will nonetheless retain certain cases 
of ^ as our solutions, even in the ring geometry, by 
assuming that the ring circumference 2L is much larger 
than the soliton size 1/k. In this limit the unavoidable 
discontinuity in dx4' &t x — xq — ±L is only of order 
0(e-2«^) 1, which we consider negligible. (It should 
therefore be noted that for fixed system size 2L, we do 
not treat cases with k arbitrarily close to zero; but we do 
treat any fixed k, for sufficiently large L.) Throughout 
this paper, whenever we discuss the case with periodic 
boundary conditions, we will present equations that are 
only strictly valid up to corrections of this order, without 
further comment. 

The grey soliton solution is continuous, and further- 
more periodic on the ring, if the phase difference across 
the whole system is Ac/i^ = 2rn7r with integer m. Choos- 
ing the chemical potential /i and the relative soliton speed 
/3 as the two independent parameters determines the al- 
lowed values of the gas velocity 



27rTO 



2L 



larctan^ 
L K 



(11) 



It follows that V has discrete values, and for finite L, the 
gas velocity cannot actually be zero for any \/3\ < c, but 
has a minimum value of order 1/L. Since the resulting 
phase gradient of is constant over most of the entire 
2L circumference, it makes a finite contribution to the 
total angular momentum. For some purposes this must 
be borne in mind, though for others a w of order 1 /L will 
be entirely negligible. We may still consider v arbitrarily 
close to any given value, for sufficiently large L. 



Inserting (|T2)) into our Hamiltonian ([6]), we find that 
terms linear in vanish because ipQ satisfies the time- 
independent GPE. The second order term H2 in H is 
therefore leading, and we ignore all higher order terms in 
this paper: 



i/2 = ^ IdxSr 



(13) 



The canonical equation of motion derived from H2 for 
6ip co-incides exactly with the linearization of the GPE 
about -ipo- 



idtStp 



1 



di ^ i{p - v)dx + 2|Voi' + i^lH* 



(14) 

This is the time-dependent Bogoliubov-de Gennes equa- 
tion (BdGE). Solving it in terms of temporally harmonic 
normal modes is equivalent to diagonalizing H2 by a 
canonical transformation. Our contribution in this pa- 
per is to do so explicitly and exactly for soliton back- 
ground ■00 for all values of /3. This expands upon the 
previously known solution for 13 — Q [23, , .24] , and signifi- 
cantly extends the set of backgrounds for which explicit 
BdGE solutions are known. In this section we set up the 
problem in detail, and discuss the orthonormality and 
completeness of BdGE normal modes, before presenting 
our explicit solutions in the next sections. 



Canonical normal modes 



We begin by expanding 



5i^{x,t) = J2iRv{x,t)qr,{t) + iS^{x,t)p^{t)) (15) 



IV. LINEARIZATION OF THE GPE 

We can obtain a larger class of solutions to the GPE 
than ([7]) alone, in the form of solutions that are close to 
it: 



ipix, t) — ^po{x) + 5'ip{x, t) 



(12) 



with (5-0 being small. In the Bose gas context, an 
operator- valued version of this perturbation 5^ is used to 
describe the leading quantum corrections to the Gross- 
Pitaevskii mean field theory. Within the linearized dy- 
namics we will discuss, however, the quantization of 
(5'0 is trivial, being that of a non-relativistic free field. 
We will therefore continue classically, analyzing the lin- 
earized Hamiltonian field theory of Siplxjt). Application 
to the quantum problem can be achieved with literally 
no greater change than simply placing operator accents 
" over all of our classical canonical variables. 



for a set of real-valued classical canonical variables 
{'ZrjjP))}, labelled with index rj. The summation over 77 in 
(|15p is used in the symbolic manner. For the finite system 
it is a true summation over discrete values of 77, whereas 
for an infinite system it stands for the summation over 
discrete values and an integration over continuous values 

of T]. 

Eqn. (IT5|) defines i?,, , S*,, , g,; , , given Sip- It is to be 
understood as a mapping from one set of infinitely many 
degrees of freedom, namely the real and imaginary parts 
of ijj for each position x, to another set, namely the pairs 
qrnPri for each r/. In order to be sure that we have fully 
solved the evolution of Sip, we must require that this 
mapping be complete; and we would additionally like 
the mapping to be one-to-one. This means that every 
possible tp can be obtained with one and only one set of 
{QqiPv}- Hence, (jlSp should be invertible, and for any 
instant t we can regard each q^j and as a functional of 
ip{x). 
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Since canonical co-ordinates are important for many 
applications (including quantization), and since any so- 
lution of a linear dynamical system can be expressed very 
easily in terms of phase space co-ordinates, we can finally 
demand that the set {qriTPr/} be canonical variables. Con- 
sequently, they fulfill the fundamental Poisson bracket 
relations 

[qn^Pn']=^W [PniPv']^iln^<ln']^0 ^ (16) 
where the Poisson bracket is defined as 

6A SB 6B 5A 



[A, B] = -ij dx 



(17) 

for any functionals A, B. (Throughout this paper the 
delta with mode indices is a Kronecker delta for discrete 
degrees of freedom and has to be understood as a Dirac 
delta distribution for variables that depend on a contin- 
uous index.) 

The definition (jl7p immediately implies that 



[6ilj{x),S^*ix')] = -i6{x-x') 
[Si;{x),6i^{x')] =0. 



(18) 



Inserting (ITSI) into (|18l) and applying ([TE)) yields the com- 
pleteness relations (or partition of unity) 



n 

Y,{Rvix)Sr,{y)-S,,{x)R^{y))=0. 



y) 



(19) 



These latter must be checked to confirm that all normal 
modes have been found, and none was overlooked. 

We show in Appendix O without any reference to the 
BdG equations, that invertibility and canonicity of (jl5p 
imply the orthonormality relations 



(20) 



J dx {R,jR*, - R*R,y) = , 



which in turn imply that the Poisson bracket is equiva- 
lently given as 



dA dB dA dB 



(21) 



as indeed it may be given in terms of any set of canon- 
ical variables. This means that orthonormality of BdG 
modes does not need to be understood as a fortunate co- 
incidence due to the particular form of the BdGE, but can 
be considered a general consequence of Hamiltonian dy- 
namics, which the BdGE of course obey. And neither the 



completeness nor the orthonormality relations are conse- 
quences of quantum mechanics; the correspondence of 
classical Poisson brackets, and quantum commutators, is 
exact, with a factor of ih. Both these relations stem sim- 
ply from Hamiltonian mechanics, which is common to 
both classical and quantum physics. It follows from the 
BdGE (gS]), though, that the inner product in 1^ au- 
tomatically vanishes, for all cases where 7^ flf^,. This 
will be addressed in Appendix [X] 



B. Time-independent Bogoliubov-de Gennes 
Equations 

We now turn to the particular form of the BdGE equa- 
tions. Anticipating that the normal modes of a linear 
system evolve as harmonic oscillators, we assume that 
the time evolution of the canonical variables is derived 
from a Hamiltonian of the form 



2^ \ m 



(22) 



in which is a mass which will be dimensionless in our 
chosen units, and is the eigenfrequency of the normal 
mode labelled by rj. The Hamiltonian evolution gener- 
ated by this form of H2 implies a set of linear equations 
among the _R^(a;), S'^(x), namely the time-independent 
BdGE: 



m^nls^ix) = HbR^{x) + ^IjIr;{x) 

■m:^^Rr,{x) = HbS,-i{x) - ■iIjoS*{x) 



with the differential operator 



Hf 



1 



dl + iip-v)d, + 2\%\-'-fl 



(23) 



(24) 



Inserting these eigensolutions i?,, , S*,, into and 
will reveal that H2 as given by is indeed equal to 
((22)) . This can be seen simply by applying (|23)) and 
the orthonormality properties, and integrating by parts, 
without needing to use the explicit forms of i?^ , S'^ that 
we supply below. This confirms that the which 
solve ([23)) are those which diagonalize H2, and that the 
solutions to the time-dependent BdGE (|T4)) are indeed of 
the form psp . with the explicit time dependence 



(25) 



g,,(<) = 9i;(0) cosfijjt + P^^^ sin ftnt 



Pri{t) = P))(0) cos Vlrji — m,ii},iq,i{0) sin 



Cases with O,, — > are covered here simply by taking the 
limit. 



C. Normal mode spectrum 



It is possible in general for Bogoliubov eigenfrequencies 
fl„ to be complex. In this case the system is dynamically 
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instable, and although it implies that the linearization 
approximation must break down on logarithmically short 
time scales, within that time frame it is by no means un- 
physical. The case of complex il,, requires several signifi- 
cant generalizations to the Bogoliubov-de Gennes theory 
for real f2^, but in the problems treated in this paper, it 
turns out that all il^ are indeed real. We prove this sim- 
ply by demonstrating that the set of modes with real 17^ 
is complete in the sense of (|19p . We can therefore discuss 
the nature of the BdG normal mode spectrum without 
treating the extra complications that arise for complex 

Since (l23l) is a coupled system of two complex sec- 
ond order differential equations, there are in principle 
eight linearly independent solutions for each eigenvalue 
flrj. The overall normalization of any solution is fixed, 
up to a phase, by the orthonormality condition (1201) . 
The phase is arbitrary, but in a slightly non-trivial way. 
Because of the explicit complex conjugates appearing in 
the equation, rescaling a solution (i?^, S'^) e*"(i?,,, 5*,,) 
does not in general provide another solution. The phase 
freedom that does exist is rather 



i?^ — > R,j cos a — im^jfljjSri sin a 
R 

Sjj S,i cos a — i p— sin a . 



(26) 



Considering ([15]) and ([22|) , it is easy to see that the second 
solution merely reproduces exactly the same 5ip as the 
first, but with the co-efficient relabeling 



Qrj ^ qri COS a H p— sm a 

Pr] — > Pr] COS a — m^il^g^ sin a . 



(27) 



This relabeling happens to be a canonical transforma- 
tion; but this does not change the fact that these two 
solutions do not represent distinct degrees of freedom in 
(5'0, but merely a double counting of the same configu- 
ration. Eliminating such double counting by picking one 
value of a arbitrarily is equivalent to discarding half of 
the eight solutions, namely those given by a-t-7r/2. There 
are thus only four distinct canonical solutions for each 
flri- Boundary conditions such as square-integrability or 
periodicity will generally eliminate at least two of these, 
providing zero, one, or two normal modes per non-zero 

A further symmetry is relating solutions with different 
niri- By canonical rescaling of q,^,Pri it is possible to set 
\mri\ — >■ 1 without loss of generality - and for the rest of 
this paper we will assume that this has already been done. 
Please note that it is not possible to change the sign of 
m,,. Both signs are possible, and both must be considered 
for any fi^, because if [R^, Srj) is a solution to (P5)) . then 
{Rff, —Srj) is a solution for — m,,. Examining ([TC|. 

however, we see that these two solutions correspond to 
exactly the same configuration of dip, with the trivial co- 
ordinate relabeling — >■ —pri- They therefore represent 
the same degree of freedom, counted twice, and we not 



only can, but must, discard one and retain only the other. 
Inspection of (PO]) shows that the co-efficients qrj,Pr^ will 
be canonical for one and only one solution in each of these 
pairs, since [q, —p] — —[q,p\ trivially. So we can maintain 
canonicity by always choosing the canonical co-efficients 
to retain. 

This means that the sign of for each mode is just 
as much a physical fact as its eigenfrequency il^. The 
sign of flrj itself has no meaning whatever; it has no ef- 
fect in either the time-independent BdGE (|23|) or the 
time-dependent solution (P5|) . fi^ may therefore be taken 
as positive without loss of generality. Inspecting H2 in 
(im , one sees that it is the sign of m,, which determines 
whether the energy contribution from exciting a given 
mode is positive or negative. A negative mass simply 
means an energetic instability; an excitation of this mode 
lowers the energy. The existence of such modes is of 
course incompatible with strict thermodynamic equilib- 
rium, but it is perfectly physical for perturbations around 
an excited background such as our grey soliton, which in 
many systems may nonetheless be very long-lived. 

We remark finally that in many BdG cases the index 77 
can immediately be identified with a Fourier wave num- 
ber k, but in a problem which is not translational invari- 
ant such as ours, this is in general not so. However, in 
our soliton problem there is an asymptotic translational 
invariance, and most of our modes will indeed have a 
Fourier index 77 — )■ fc, such that for k\x — xo\ 3> 1 the be- 
have as sinusoidal functions of kx. But as will be shown 
later, there will be one extra mode of frequency zero. For 
this mode we will use the single additional index rj = z. 



D. Comparison with complex co-ordinates 

In the literature on the dilute Bose gas in particular it 
is more common to employ a different form of expansion 
from (|15p. which is, however, fully equivalent to it. This 
simply uses complex variables a,, , a* instead of the real 
pairs Qn^Pn- 

Si; = J2iu^{x, t)a^(t) + v;{x, t)a*^{t)) . (28) 



For each individual non-zero frequency the transforma- 
tion between the sets is defined by 











i 

+ —= 


Pn 


) 


a; 


V2 ( 






-P,i 





(29) 



The two variable sets are connected by a so-called ex- 
tended canonical transformation [25} with a scaling pa- 
rameter X — ~i. In the quantum version of this prob- 
lem, just as the "0 field becomes a destruction operator, 
and our phase space co-ordinates qrj,P7j become operators 
QrnPr] with canonical commutation relations, these com- 
plex co-ordinates a^, a* become annihilation and creation 
operators a^, aj^ for quasi-particles. In this sense there is 
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nothing particularly quantum mechanical about the Bo- 
goliubov expansion ([25)) . any more than there is about 
(fTS]) . They are simply equivalent alternative representa- 
tions of linear Hamiltonian dynamics. 

Expressing the non-zero frequency part of (j22p in terms 
of the complex variables yields the familiar Hamiltonian 



H2 = Y.n„a;it)a^{t) 



(30) 



In analogy to ([2^ . the resulting linear equations in the 
complex co-ordinate representation are 



(31) 



with the same differential operator Hb defined in ([24)1 . 
This is a form of the BdGE nowadays widely used in 
the literature. In the a, a* formalism it is conventional 
to make the sign of fi,, meaningful and eliminate m,j by 
absorbing m,( = ±1 in f2^. 

In this paper we will mainly use the representation 
with the set {qr,,Prj}, because although it is in some re- 
spects more cumbersome, it has one important advan- 
tage: since the Hamiltonian ()22p does not vanish for zero 
frequency, unless happens also to be infinite, there 
is no ambiguity or singularity in the {(?,,, Pr;} expansion 
when the mode frequency SI,, — > 0. In the {a^,a*} rep- 
resentation, in contrast, a mode for which fi,, vanishes 
must either have its Hamiltonian exactly zero, or else 
have some singularity in (j28|) . To avoid such complica- 
tions we retain the {qrijPri} representation. We will be 
devoting considerable discussion to the modes with zero 
frequency, both because their technical subtleties require 
it, and because one of them is of particular interest, in 
that it represents the motion of the soliton itself. 

We are now ready to present the main result of this 
paper: the explicit solution of ([23]) for a complete set of 



V. SOLUTIONS OF THE BOGOLIUBOV-DE 
GENNES EQUATION 

A. Non-zero frequency solutions 

In this subsection we present the explicit non-zero fre- 
quency solutions of the BdGE (for which 77 — > A;) for the 
general case of a grey soliton; their derivation through a 
factorization technique analogous to Schrodinger super- 
symmetry is shown in the next subsection. The special 
case of the zero frequency solutions follows further below. 

If we rescale all the positive masses to mk — 1, the 
solutions to the Bogoliubov equations (|23p for non-zero 
frequency in the co-moving frame in case of both finite 



and infinite system can be expressed as 

Rk{x) = e-'"^ ^/rmRerkix) + ilmskiS;)) (32) 



Sk{x) 



(Resfc(i) + ilmrk{x)) , 



(33) 



where we have introduced the abbreviation x = x — xq. 
The functions rfc, are given by 

rk{x) = -Jie*'^"^ — - 2f3nk+ik^Kt&nhKX 
"fc V 2 

-I- kn'^sech^Kx^ (34) 
Sfc(x) = Nfc e*''^(fc + i2KtanhKx) (35) 
with the mode frequencies flk being 



ftk + f3k 



k^ 

T 



c?k'^ 



(36) 



This represents exactly the same dispersion relation as 
for a constant background solution = e~™^(i/3 ± k) 
equal to the grey soliton background solution in the limit 
X — >■ ±cx). For a dark soliton (/3 = 0) this has already 
been demonstrated in [H, . 

The normalization constant that ensures ([20)) is 



k 



(37) 



for delta-function normalization in the infinite domain, 
with a continuous spectrum of fc, and 



=*^8L(17,+/3fc)-4^ 



(38) 



for the finite ring with its discrete spectrum. For the 
demonstration of the orthonormality of these modes, see 
Appendix [X] 

For the finite ring geometry it is also necessary to se- 
lect the discrete set of k for which the solutions are peri- 
odic (up to exponentially small corrections, as discussed 
above). Just as with the background solution i/jq, these 
are obtained by imposing continuity for the 5^ modes 
and their first derivatives. As long as ip^ is thus periodic, 
the condition for periodicity of 5^ reduces to 



k cot kL — 2 



nk 



(39) 



In case of the dark soliton /3 = 0, this condition is equiv- 
alent to the one given in [24]. 

From the explicit solutions above one can see that any 
perturbation of the grey soliton background made from 
the normal mode spectrum passes through the soliton 
without being reflected, which is expected for the intc- 
grable model [ll|3l|. 

The proof that the above non-zero-frequency modes 
together with the zero frequency modes described in Sec- 
tion IVII below are a complete set in the sense of (|19l) is 
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given in Appendix |B] It proceeds simply by direct in- 
tegration, which can be performed analytically, but is 
somewhat involved. It is not hard to see that no ad- 
ditional finite-frequency modes are expected, however. 
Due to the symmetry of the Hamiltonian for non-zero 
frequency modes mentioned in IIV( it is only necessary 
to find four independent solutions for a given eigenfre- 
quency. For non-zero frequency fl, the Bogoliubov dis- 
persion relation (j36[) is a quartic equation in the mode in- 
dex k. Assuming all frequencies are real and that ^ 0, 
this equation must in general have two real roots and two 
mutually complex conjugate roots for k [s^ . Since the 
complex k solutions are un-normalizable in the infinite 
domain, and can be shown incompatible with periodicity 
on the ring, we are left with two distinct solutions with 
real k for each ft, and these are the solutions we have 
displayed. 



B. Comparison with constant background solutions 

It is quite instructive to redo the complete linearization 
procedure for the constant background solution = 
e~™^(i/3 ± k), which is equal to the soliton background 
solution ^/jQ at a; — >■ oo or X — >■ —oo. As mentioned before, 
one obtains the same dispersion relation p6l) , but instead 
of ([M)) and the solutions are given by 



ikx 



(40) 
(41) 



The normalization constant of these non-zero frequency 
solutions is the same for the infinite and finite system 
and is given by (l37l) . Except for the phase factors the 
solutions for the BdGE for constant background are con- 
stants. Comparing these with (p4l) and (|35|) one immedi- 
ately sees that solutions for the constant background 'ip± 
coincide with the soliton solutions at J — > ±cxi. 



rewrite the BdGE ([23)1 in terms of the functions , Sfc , 
related to the Rk , Sk by: 



Rk + ^kSk 
Rk — flkSk 

With the differential operators 



rk 



Sk 

si 



1 



2Ktanh kx) 



= —j=[—dx + 2KtanhKi) 



(42) 



(43) 



and the explicit form of the grey soliton background, the 
BdG equations (|23| can be cast into the alternative form 



Vlk 



rk 




Sk_ 





iV2f3Q QQ^ 



f 2(/32 - ^2y 





rk 




Sk 


X 





(44) 



The displayed functions depend on x = x — xq. If we 
define a new function E by Q^S = Sk, one obtains from 
the bottom line of (l44t 



{nk+iV2l3Q^)^ = Q 



rk 



(45) 



Acting with Q from the left on the top line of (|44)) and 
using the relation (|45l) and definitions P5|) . the top line 
finally yields 



i2i3nkd. 



4^ 



«'a2)I] = 



(46) 



Although this fourth order equation does not contain hy- 
perbolic functions and looks very simple, it is equivalent 
to the coupled set of second order equations (|44)) ; we have 
not made any approximations. ()46|) has the obvious solu- 
tion Yik — Ce'*^^ with some normalization constant and 
the Bogoliubov frequency fife equal to (|36p . Applying the 
differential operator to directly gives us Sk, and 
lets us obtain by solving the first order equation (|45p . 
This verifies the non-zero frequency solutions p4[) and 
(1551) of the preceding subsection. 



VI. ZERO FREQUENCY SOLUTIONS OF THE 
BDGE 



C. Derivation of non-zero frequency solutions via 
factorization 

In one dimensional linear Schrodinger systems, scat- 
tering off a potential can sometimes be simplified with a 
factorization of the Hamiltonian. Factorization is a tech- 
nique which relates a problem at hand to a system with 
a simpler potential by using a proper superpotential. In 
this way one can easily solve the one dimensional scatter- 
ing due to a potential of the form sech^ kx. for example, 
since it is connected to a constant potential [331 ■ have 
demonstrated in the preceding subsection that the grey 
soliton background acts like a reflectionless potential for 
the non-zero frequency modes just like the K^sech^KX- 
potential, which suggests the usage of a method similar 
to the factorization method. To do so, we are going to 



We now come to the zero frequency solutions, often re- 
ferred to briefly as zero modes, which present special dif- 
ficulties but are particularly important from some points 
of view. The zero frequency solutions solve the coupled 
system of equations 

- HbR^ + i'lR; (47) 
m-ii?„ - HbS,j - i>ls; (48) 

We will begin with the solutions for the infinite system 
without boundary conditions, and then derive the solu- 
tions for periodic boundary conditions. We will end this 
Section by discussing an interesting invariance property 
of the zero modes. 
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A. Infinite system witiiout boundary conditions frequency modes (see Appendix [P]) : 



A straightforward way to obtain solutions for the cou- 
pled set of equations (H71) and is to take the deriva- 
tive of the GPE with respect to the five parameters 
{l3,v,xo, 11,9) of the background solution. All GPE 
solutions within this five-parameter family are time- 
independent, up to a uniform phase prefactor e"*''*, 
and in some co-moving frame. Small perturbations 
within this family will therefore automatically have zero 
frequency. Those perturbations which do not shift /i 
or the lab-frame soliton velocity will be strictly time- 
independent. Those which do infinitesimally perturb 
fl will produce linear time dependence dip oc t, since 
g-zSp.t ^ 2^ _ iSp,t. Perturbations of the soliton's speed 
similarly induce a linearly time-dependent translation of 
the soliton, in any fixed frame. And linear time depen- 
dence is also zero frequency. 

These variations do not provide all zero-frequency 
BdGE solutions, but they do provide all the solutions 
needed for completeness. We can confirm this by finding 
the remaining linearly independent solutions, by other 
methods, and showing directly that they are unphysi- 
cal. It is important to remember that for the infinite 
domain all five parameters are independent, apart from 
the condition k > 0; the periodic boundary condition 
that we subsequently consider will select a subspace out 
of this solution space, reducing the number of physical 
zero modes. 

We demonstrate the procedure with an example. In- 
serting the background solution ^ into the GPE (jH) and 
taking the derivative with respect to 9, yields 



1 



(49) 

Hence, we have found a first solution to the homogeneous 
equation (|T7)) : 



Ri = iipo ■ 



(50) 



Taking the derivative with respect to the other solution 
parameters xq,P,v and /i, keeping in mind that we have 
abbreviated p, = fi + v(3 — does not lead to all four 

solutions to ()47|) . however. Varying jl oi (3 — v produces 
(linearly) time-dependent perturbations, unless we also 
change frame and the e~*'^* prefactor in ^ . Hence these 
perturbations do not yield further solutions to (j47p . Per- 
turbations which do not affect /i or /3 — w, though, give 
two more solutions to the homogeneous equation: 

i?2 = {iv - dxo)ipQ = {iv + dx)tpo (51) 
^3 = -dpipo - dvi}o + l^d^ipo . (52) 

Since (|Tfl) is a system of two coupled second order 
linear differential equations, for the real and imaginary 
parts of R, there must exist a fourth linearly independent 
solution. We have obtained this by a variant of the fac- 
torization method used in subsection IV CI for the finite 



i?4 = e ™^(3k5 sech^Ka; + StanhKx + sinh2Ki 



(53) 



with X — X — xq. 

Turning now to the equation for the zero- modes' Sri, 
we can first identify four homogeneous solutions, which 
make the right-hand side of pS)) vanish. All four are 
immediately apparent: S — >■ iRj. 

Then, for each of the four Rj on the left side of (|48p . 
there must exist corresponding particular solutions Sj to 
the inhomogeneous equation. Here we can make use of 
perturbations of fl and jS — v, which generate Sip oc t. 
Inserting the general zero mode time evolution 



Pvit) =P»,(0) 



(54) 



into Sip = R-qqri + iS'jjP,,, we can observe that perturba- 
tions to jl and l3 — V must correspond to some combina- 
tions of zero- mode Pr;(0). It follows that derivatives of V'o 
with respect to /l or /3 — w will provide the components 
Sjj — —idStp/dpn of the zero mode {R,S). We are also 
free to include any combinations of the homogeneous so- 
lutions iRj. Two simply expressed particular solutions 
that match Ri and R2 are respectively 



Si 
S2 



idpipo 



(55) 
(56) 



It can then simply be checked that {Ri, Si) and (i?2, ^2) 
are indeed zero frequency solutions of the coupled BdGE 

dUl) and dm). 

There are another two particular solutions to (|48l) . 
with i?3 and R4 on the left-hand side. In the appendix IdI 
we derive the asymptotic behavior for x — ±00 of these 
solutions 5*3 and S4, which is (as we will see) all we need 
for their consideration. 

In this way we have found all eight linearly indepen- 
dent zero mode solutions, at least asymptotically. They 
are of the form {R, S) {Rj,Sj) and {R, S) (0, iRj), 
for j = 1,2,3,4. We can exclude most of these eight 
solutions as Hamiltonian modes, however. {Rj,Sj) for 
j = 3,4 and (0,zi?4) are all unphysical, because one can 
prove from their asymptotic form that no linear combina- 
tion of them can be normalizable according to (PU)) (they 
diverge exponentially at infinity). And we may discard 
without loss of generality the two zero frequency solu- 
tions (0, iRj) for j = 1, 2, because their effect in Sip{x, t) 
could be absorbed by a shift of the canonical position 
variables qi and q2. So, in case of the infinite system 
without boundary conditions, we are left with three pos- 
sibly physical zero frequency solutions to the BdGE. 

Two of these remaining three solutions may be ex- 
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pressed as the linear combinations 



5o± 



1 



(27r7±c) 2 



7± 



Ri 
Si 



± 



R2 
S2 



± 





iRs 





1 


'-R2 




~ 2^^ 


^2 



(57) 

Here 7± = c =F /? is the hmit of r2fc/|fc| for k — )■ 0^. 
One can show that ([57]) defines two positive mass zero 
frequency solutions which are (5-normalizable. These two 
modes are precisely the limits /c — >■ 0^ of the continuum 
of non-zero frequency solutions p4l) and psp . Therefore 
they do not need to be separately added or considered; 
they are already fully represented within the continuum. 

In addition to these two positive mass modes we have 
a third linear independent mode, which we can take to 
be (i?2, 5*2 )• This has finite negative norm, and so we can 
identify 



(58) 



as a positive norm discrete mode with negative mass. 
This is the soliton translation mode, as one may see by 
recognizing that its effect in Sip is precisely to make a 
small translation of the soliton. 

We prove in Appendix [B] that the continuum solu- 
tions alone are not complete, but that adding {Rz,Sz) 
as one additional discrete mode provides completeness. 
We can therefore be confident that we have obtained all 
modes of the system as BdG solutions. There remains 
one point, though, which may seem to call for comment. 
Although we have identified three linearly independent 
zero frequency solutions, the three are not all orthogonal 
to each other according to our inner product as intro- 
duced in ([^ . One source of this curious difficulty is 
the fact that {R, S) — > (0,ii?3) has zero norm according 
to the inner product. The reason that this subtle issue 
raises no problems for completeness is simply that the 
zero modes contained within the continuum are a set of 
zero measure. For any finite frequency, however small, 
there are two normalizable modes, and not three; and so 
the 'extra mode' (0,zi?3) at exactly zero frequency has 
no effect. 

We attribute the appearance of the extra zero-norm 
solution to the fact that the boundary condition, which is 
only that R and S should not diverge too fast at infinity, 
is somehow too weak. If we replace the strictly infinite 
line with a periodic ring, however large, then the (0, iRs) 
solution ceases to be independent of the other two, but 
instead must be added to them in particular proportions, 
in order to meet the periodic boundary condition. 



B. Periodic boundary conditions 

Before deriving the periodic zero frequency solutions 
by taking derivatives of the GPE as in the preceding 
subsection, we want to emphasize two details. Firstly, we 
remind the reader that the results for the periodic system 
presented in the following are valid to order C'(e~^''^), 



which is an excellent approximation for kL ^ 1. Sec- 
ondly, the gas velocity v, the relative soliton speed /3 and 
the chemical potential /i are not independent parameters 
for the periodic system, but must together satisfy a pe- 
riodicity condition, as explained in Section IIIII It turns 
out that cumbsersome expressions simplify most easily if 
we choose the relative soliton speed /3, and the chemical 
potential fi to be the two independent variables, in terms 
of which V may be expressed. Hence, derivatives with re- 
spect to /3 and fi are now different from the infinite case, 
because v must now also be varied whenever /? and fi are. 

Inserting the grey soliton background ([7]) into the GPE 
(HI), and differentiating this equation with respect to 
the remaining four independent parameters (6*, xq, /3, /i) 
yields the four constituent functions i?i,2, 5'i_2 of two so- 
lutions {R, S) — >■ (.Ri,2, 'S'i,2)- Since the background solu- 
tions within this four parameter family are all periodic, 
all their derivatives with respect to these four parame- 
ters are automatically periodic as well. A linearly inde- 
pendent basis for these two periodic zero modes can be 
expressed as two different periodic linear combinations of 
the three non-periodic solutions for the infinite domain 
considered in the previous subsection IVI Al 



'Ri' 

Si. 




Ri 




'R2 

S2. 




'R2 

.^2. 





i?2 - vRi 
S2 - l3Si + zi?3 

i?2 - vRi 
S2 - /35i + ii?3 



(59) 
(60) 



The explicit expressions for d^v and dpv can be found 
by differentiating (fTTj) with respect to the corresponding 
parameters and are both of order 0{1/{kL)): 



2KLfi 



(61) 



Besides these two periodic zero frequency solutions 
(i?i_2, 'S'1,2), and the corresponding homogeneous solu- 
tions (0,«-Ri,2) that can be absorbed into shifts of qi^2, 
there must exist four more zero frequency solutions. 
Since in comparison with subsection IVI Al above we have 
changed only the boundary conditions and not the equa- 
tion, these solutions must simply be (i?3,4, 5*3,4) and 
(0,zi?3,4), where {Rj,Sj) are given in Appendix [Dl It 
is not difficult to show that no combinations of these 
four solutions can be made smoothly periodic, and so 
(i?i^27 »S'i,2) span all the zero modes that we need to 
canonically represent Stpix^t). 

It remains only to select combinations of these solu- 
tions, and possibly perform the mass-changing transfor- 
mation {Rj, Sj) — ^ {~Rji Sj), to assemble from our two- 
solution basis a pair of positive norm solutions that are 
orthogonal to each other. This is accomplished with 



i?0 


i 


_Sq 




'Rz 


1 


Sz_ 


~ 2v^ 



Ri 
Si 

-R2 
S2 



o 



1 



o 



(Ki)3/2 
1 

kL 



(62) 
(63) 
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Here the corrections of higher order in 1/{kL) may be 
worked out exphcitly, but we will assume that L is large 
enough to ignore them. The negative mass solution (|63|) 
is equal to the corresponding solution ([55]) for the in- 
finite system, up to the indicated order, so we use the 
same subscript z. It is the zero mode associated with 
the translation of the soliton. The periodic positive mass 
solution (|62p . labelled with subscript 0, has been normal- 
ized and can be identified with the usual phase transla- 
tion zero mode, which is present for any stationary Gross- 
Pitaevskii solution, to leading order in 1/{kL). 

C. Symmetry of the zero modes 

The main result of this paper has now been achieved: 
we have presented a complete set of BdG solutions, in- 
cluding zero modes, for any grey soliton background. The 
only remaining question is whether our solution is unique. 
Since all our non-zero frequency modes are in degenerate 
pairs, it is obvious that arbitrary linear combinations of 
each degenerate pair may be taken; this trivial freedom 
is the only one remaining for the finite frequencies. 

The two zero modes of the finite system, namely the 
soliton translation and the phase translation mode, can 
also be mixed into each other, but the transformation 
is a bit different. It is straightforward to show that the 
BdGE ([47)1 and (|48)l are also solved by the following new 
zero-frequency solutions for any real 9: 







cosh 6' 


— sinh 9 





" 




'Rq' 






— sinh 9 


cosh 9 










Rz 


S'o 










cosh6' 


sinh 9 




So 


Is'J 










sinh 9 


cosh 9_ 




Sz. 



(64) 



The two solutions are again of opposite mass. Consid- 
ering and ([22]) one can see that this transforma- 
tion leaves Sip unchanged, and is equivalent to a proper 
canonical transformation of the (71,2, Pi. 2 for the two zero 
frequency modes. 

While this symmetry is thus in one sense trivial, it is 
also surprising. It mixes the soliton motion and global 
phase translation degrees of freedom; and moreover these 
are modes of opposite mass. One may be tempted to sup- 
pose that positive and negative mass modes could never 
be confused, because they would have drastically differ- 
ent behavior. But in fact they may not always be so 
distinct, after all: the zero modes can be mixed by the 
canonical transformation associated with (j64[) . 

VII. DISCUSSION 

This has been a technical paper on a technical subject. 
Since there are not many exact but non-trivial BdG so- 
lutions known, we have tried to use our solution to pro- 
vide a pedagogical discussion of important aspects of the 
general BdG problem. Many of the properties of our 
solutions, including the fact that we were able to find 



them, seem to depend on particular detailed properties 
of certain hyperbolic trigonometric functions, and in this 
sense are purely technical features of a special case. It 
is in particular disappointing to report that, although 
the grey soliton GPE solution itself may be extended 
into a class of multi-soliton solutions given exactly by 
certain elliptic functions, and this class may be shown 
to include all time-independent GPE solutions with uni- 
form external potential in one dimension, we have been 
unable to extend our exact BdG solutions to these cases. 
The solutions we have found seem to depend crucially on 
properties of tanh and sech^ that are not shared by their 
corresponding elliptical generalizations. But the problem 
of BdG normal modes in a grey soliton background is in- 
teresting from some rather fundamental viewpoints, and 
we will close our paper by discussing two of them. 

A feature of the stationary grey soliton solutions which 
has on occasion caused excitement is the fact that they 
seem to include sonic event horizons. The local speed of 
sound in the hydrodynamic approximation to the GPE 
is 1-01; the local fluid speed is dxSc^gi^- It is easy to 
show that this local fluid speed exceeds this local sound 
speed within a finite range — a;o| < Xh for a certain j3- 
dependent Xh , apparently providing a black -I- white hole 
scenario with horizons at a; — a;o = ±a^/i . The very prob- 
lem we have solved would therefore seem to be an ideal 
candidate for terrestrial investigations of the Hawking ef- 
fect in black holes, by examining quasiparticle pair pro- 
duction after quantization. Unfortunately, however, our 
exact results show that this does not work: grey solitons 
do not really have horizons. No perturbation wavepack- 
ets are reflected at any point in the soliton background, 
regardless of which direction they move in; modes of any 
wavelength can freely pass through the supposed hori- 
zons, in any direction. 

This disappointingly ordinary behavior illustrates an 
important caveat which must never be overlooked in con- 
structing black hole analogs in fluids. Sound waves only 
imitate light, with a maximum speed and no dispersion, 
in the hydrodynamic regime. This means that the wave- 
length must not be too short, so that the hydrodynamic 
picture does not break down. Yet even for real black 
holes and real light, the horizon concept is only applica- 
ble for wavelengths short in comparison with the length 
scale over which the background changes, so that the 
geometric optics limit in which light propagates on null 
geodesic rays becomes valid. Light with a wavelength 
much longer than a real black hole's Schwarzschild ra- 
dius is not trapped by the event horizon, but simply 
diffracted. The requirement for a sonic horizon is there- 
fore that there exist a range of wavelengths which are si- 
multaneously long enough to be hydrodynamic, but short 
enough to be geometric. For a grey soliton it turns out 
that there are no such wavelengths: perturbations are all 
either non-hydrodynamic or non-geometric. Our exact 
solutions then merely confirm the fears of horizonless be- 
havior which will be raised by careful estimates of validity 
regimes. 
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Another fundamental question which grey sohton BdG 
solutions can address is the emergence of collective co- 
ordinates in classical fields and quantum many-body sys- 
tems [sS]. Such collective co-ordinates represent many 
fundamentally interesting and even practically impor- 
tant degrees of freedom in the real world. The men- 
tal degrees of freedom of conscious beings, for example, 
are presumed by substance monists to be collective co- 
ordinates of the large and intricate quantum many-body 
systems known as brains. Such co-ordinates stand out in 
some way from the quasi-continuum of other degrees of 
freedom, but they do not entirely decouple from them; 
the closest we can accurately come to ignoring all those 
less interesting other modes is to treat them as an envi- 
ronment or reservoir, and describe the evolution of the 
collective co-ordinates as a kind of generalized Brownian 
motion. This is understood to introduce thermodynam- 
ics and quantum decoherence, but understanding exactly 
how this happens and what it means is a goal still un- 
achieved. 

Perhaps the most basic question is simply, What is 
it that makes the collective co-ordinates stand out from 
the quasi-continuum crowd? One might call this very 
basic issue 'the Hamiltonian gavagai problem', after 
the philosopher W.V.O. Quine's identification of a sim- 
ilar problem in linguistics [35]. Quine imagined that 
a linguist learning an unknown language from a native 
speaker might see a rabbit run across a field, and hear 
the native exclaim, "Gavagai!" While the hypothesis that 
'gavagai' was simply their language's word for 'rabbit' 
might seem obvious, Quine pointed out that it might just 
as well mean some part of a rabbit, or some aspect of the 
rabbit's motion, or some relationship between the rabbit 
and the field, or any number of other concepts. While 
Quine's 'gavagai problem' was formulated as a problem 
in the theory of translation between languages, it has 
become a touchstone in the philosophy of meaning. A 
much simplified but basically similar problem presents 
itself in the Hamiltonian dynamics of many degrees of 
freedom. Infinitely many canonical coordinate systems 
are possible; some are evidently better than others; but 
we lack an explicit theory to determine which ones might 
be best. In particular, although it is a familiar fact that 
a small number of degrees of freedom often stand out as 
prominent, we lack a theory of prominence. 

While it is possible that such problems simply lie out- 
side the subject of physics, and will remain grist for phi- 
losophy, they do seem difficult to avoid when discussing 
important physical processes, such as dissipation and de- 
coherence, at a fundamental rather than phenomenolog- 
ical level. Perhaps it is worth looking for a dynamical 
theory of prominence which might, like adiabatic theory 
and in the study of dynamical chaos, prove both rigorous 
and rich, despite being based on concepts that were not 
originally seen as part of mechanics. 

Such a goal surely lies well beyond the kind of lin- 
earized analysis we have presented in this paper, since 
linear modes can always be exactly decoupled. Within 



linear theory, there is no obvious reason why a collec- 
tive co-ordinate should not simply be one of the normal 
modes, and hence decouple exactly from all the other 
modes, after all. Yet the linear regime may still offer 
some insights and hints. 

In the first place we confirm that the problem is subtle. 
The motion of the soliton, which is the obvious collective 
degree of freedom in this problem, is indeed represented 
as one of the BdG normal modes, namely the negative 
mass zero mode. But we have seen that there is no ab- 
solute distinction between this mode and the other zero 
mode, which is a long wavelength background phase gra- 
dient mode that is present even in the absence of the 
soliton. A canonical transformation can mix the two 
modes, while preserving exactly the form of their lin- 
earized Hamiltonian. One can conjecture that nonlinear 
dynamics, which in the GPE is local in space, may some- 
how prefer that canonical representation of the two zero 
modes which makes the soliton motional mode as spa- 
tially localized as possible. Since the soliton's phase pat- 
tern is very nonlocal, however, even though its modulus 
deformation is localized exponentially well, the most lo- 
calized zero mode possible is still not perfectly localized, 
but does include long range perturbations of the phase of 
-0. Yet since those perturbations are spatially uniform ex- 
cept near the soliton, and a uniform phase shift is in most 
contexts unobservable, it is not clear to what extent this 
represents a 'genuine' long range extension, rather than 
something like a topologically nontrivial gauge configu- 
ration. 

It is instructive to contrast the subtlety of these issues 
for the grey soliton background with their triviality in 
the bright soliton case, where the sign of the term in 
([T]) is negative rather than positive. A bright soliton is a 
localized blob with ^ ~ sech k{x — xo)- Rather than mov- 
ing through an asymptotically uniform background field, 
it is itself the entire field; far from xq , the field vanishes 
exponentially. The bright soliton can move, but its mo- 
tion is identical to the motion of the entire configuration; 
there is no possibility of relative motion between the soli- 
ton and the background, as there is the grey soliton case. 
Dynamically, this difference shows up in the fact that 
the grey soliton's motional zero mode lies at the bottom 
of a gapless continuous spectrum of other modes, while 
the bright soliton's zero mode is separated in frequency 
from the continuum by a finite gap of width (See 
Appendix lEl below.) Distinguishing soliton from back- 
ground, as a collective degree of freedom, can therefore be 
done easily and unambiguously for the bright soliton, be- 
cause of this large separation in time scales. For the grey 
soliton there is no time scale separation, but only a con- 
trast in the pairing of time and length scales. The grey 
sohton moves much more slowly than any other pertur- 
bations with wavelengths on the order of its size, and it is 
much smaller in extent than any other perturbations that 
evolve as slowly as it moves. There is therefore reason to 
hope that the role of locality in nonlinear GPE dynamics 
will indeed single out a unique soliton collective degree 
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of freedom; but there is no obvious identification based 
on time scales within the hnear regime. 

The technical difficulties in going beyond the linear 
approximations used in this paper are considerable, but 
they are joined by conceptual difficulties as well, particu- 
larly in the quantum problem. The difficult issue is back- 
reaction: since the continuum modes do interact with col- 
lective degrees of freedom, beyond the linear regime, we 
must expect that the collective co-ordinate evolution is in 
some way 'dressed' or perturbed by the other modes. Yet 
if these modes are quantum mechanical, or even if they 
are classical but subject to thermal fluctuations, then 
one seems to be faced with the challenge of somehow 
incorporating thermal or quantum fluctuations into the 
collective co-ordinate itself. This may be feasible, though 
it might imply large fluctuations or even a 'Schrodinger's 
Cat' state; but it is not clear whether it is really correct, 
or whether perhaps instead the right definition of the col- 
lective degree of freedom must in some way average over 
these fluctuations, in the same way that a temperature 
or a center of mass involves an averaging by definition. 
This issue is not purely philosophical, but may in prin- 
ciple have clearly observable consequences. In a series of 
single observations of an identically prepared quantum or 
thermal system containing a soliton, do we expect to see 
in each run a smeared out soliton in the same location? 
Or do we expect to see a sharply resolved soliton, whose 
location fiuctuates randomly from run to run? 

Nonlinear perturbation dynamics in grey soliton back- 
grounds may thus provide a simple model for issues of 
fundamental importance. By providing a complete and 
explicit analytic solution of the linearized problem in this 
case, we have laid a foundation on which that further 
study can be built. 



On the one hand we have directly evaluated these inte- 
grals, from our explicit expressions for and could 
prove the orthonormality relations ([20)1 for continuous 
index k and the discrete zero mode with index z in a 
lengthy but straightforward calculation. It is much eas- 
ier, however, to prove orthonormality for modes with 
il^ ri| using the BdGE and integration by parts. 
We demonstrate this for (IA3I): 



i{R,,\iR^] 



= m^ jdx [RniHeS^ - V'g^D* 
= jdx [SliHBRr, + V'o^;) 

-s^{HBR^ + ^PlR;y 



(A5) 



In the first line we have used the BdGE ([23]) for R^ and in 
the second line we have integrated by parts. Once again 
using the BdGE yields 



i{Rrf\iR^) = ■m^mjfQ'^{S,f\iS^ 



(A6) 



As a first result we find that for non-zero frequency modes 
the relations (|A3|) and (jA4l) are not independent of each 
other. Redoing the procedure (jA5[) for (jA3p . but now 
using the BdGE for i?^ instead, reveals 



i{R,-i\iR^) — m^m,in'^{Sri\iS^) 



(A7) 



Hence, for ri| ^ fi^^, one can deduce (|A3|) . With this 
result also (|A4[) is proven by (|A6|) . 

Along similar lines one can confirm (jA2p . In a first 
step it is easy to show that 



Appendix A: Orthonormality 

We will first show orthonormality for the infinite sys- 
tem; the finite periodic case is straightforwardly similar. 
If we define for complex functions /(a;), g{x) 



{f\g)^ dx{fg* + rg) 



(Al) 



with integration over the infinite or finite domain, respec- 
tively, then the orthonormality relations ((20)) we have to 
confirm are equivalent to 



{Rn\Si)^S, 
{R^\tR^)^0 



(A2) 
(A3) 
(A4) 



The Kronecker delta is to be understood as a Dirac delta 
function 6{k~k') whenever 77 — >■ is continuous (in which 
case all three integrals are distributions, rather than true 
functions of k,k'). Where the discrete mode is involved, 
we intend Szz — 1, and 5kz = Szk = (also as a distribu- 
tion in k). 



{R,\S^) = ^{S^\R,) 



(A8) 



It follows that the orthonormality relation (IA2p is sym- 
metric under switching 77 -s-> ^ up to the constant factor. 
In analogy to the above proof it is a straightforward ex- 
ercise to confirm (IA2[) . 

The results found for the infinite system can be trans- 
ferred to the finite system with periodic boundary con- 
ditions in the following way. The normalization constant 
Nfc has to be changed (see ([57)) ) and the delta-distribution 
in (l^ni) turned into a Kronecker-delta, while the allowed 
non-zero values of the index k are given by the quantiza- 
tion condition ()39p . The two zero frequency modes pre- 
sented in subsection lVI BI can easily be made orthonormal 
in a Gram-Schmidt procedure, and they are naturally or- 
thogonal to all modes with non-zero frequency. 



Appendix B: Completeness 

For the infinite system without boundary conditions, 
the completeness relations ()19p for the set of continuous 
Bogoliubov modes ([5^ . (1551) and the discrete zero mode 
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(j58p can be proven directly, by straightforward contour 
integrations over the continuous index k. We express the 
continuous Bogohubov modes in terms of the functions 
rk and Sk given in and (1551) . Both relations then 
involve terms of the form 



dk[rk{x)sl{y) + rl.{x)sk{y)] = 5{x - y) ~ y sech kx 



dk [sk{x)sl{y) - sl{x)sk{y)] = 
dk [rk{x)rl{y) - rl{x)rk{y)] 



(tanh Kx sech^K?/ — tanh ny sech^ra) 



+ k{x — y) sech^/tysech^Kcc . 



(Bl) 



Simplifying the results of the contour integrations to re- 
veal the above simple expressions requires some judicious 
application of identities among hyperbolic trigonometric 
functions, which the authors have successfully performed, 
but the exercise is involved and we omit the details here. 
Numerical plotting of the unsimplified results will quickly 
confirm (IBip . 

The fact that the right-hand sides of (jBl[) is not of 
the simple form required by canonicity and completeness 
proves that the continuous BdG modes are not complete 
by themselves. But it easy to see that this defect is reme- 
died precisely by adding one more mode, namely the dis- 
crete negative mass zero mode associated with soliton 
motion, (i?z, Sz)- Thus we have proven that the discrete 
negative mass zero mode and the continuous modes form 
a complete set of functions. 

For the finite system with boundary conditions the in- 
tegral over the modes becomes a sum over the allowed 
indices (cf. Eq. ((23)) plus the periodic negative and pos- 
itive mass zero mode we have found in IVIBI A direct 
proof that the summation gives a partition of unity is not 
accessible due to the complicated quantization condition 
(15^1) . But it seems plausible that we can decompose any 
function that satisfies periodicity on the ring in terms of 
the periodic functions of our complete set. 



Appendix C: Orthonormality from completeness 

Here we show that the orthonormality relations ([20)1 
are a consequence of completeness and canonicity (1191) . 
together with uniqueness, i.e. invertibility of the mapping 
from {(jr;,P?j} to 5'ilj{x). 

We begin by defining the generalized matrix 



Rr,{x) 

r;{x) 



iSr,{x) 

~iS*{x) 



(CI) 



Here 77 is a Bogoliubov mode index; summing over 77 is to 
mean summing over all discrete Bogoliubov modes fc, z 
in the periodic case, or integrating over continuous k and 



then adding the discrete zero mode ry = z in the infinite 
version of the problem. We call Ax^^i a generalized ma- 
trix simply because we want to consider x and rj as two 
matrix indices, even though x, and possibly 77 as well, is 
continuous. In a similar sense we define also the general- 
ized matrix 



B„ 



[x) Srj {x) 

-iR^ix) iR;{x) 



(C2) 



These two matrices A and B are constructed so that the 
'matrix product' AB gives the completeness relation: 



'Six ~y) ■ 
5{x - y) 



(C3) 



where as usual we mean the sum over 77 to be interpreted 
as an integral in case 77 is continuous. 

Orthonormality, on the other hand, concerns the ma- 
trix product BA: 



dx Bfj^xAx,^ 



(C4) 



Simply working out the matrix product here shows that 
the Bogoliubov orthonormality property, which we are 
trying to demonstrate, is just that M is in fact the iden- 
tity matrix. This we now prove, in two steps. 

First we show that M is a projection operator, and so 
all of its eigenvalues must be either or 1. To see this 
we simply compute 

X 

= dxdy"^ Brj^xAx^xBx.yAy^^ 
= Jdx dy B,^^xNx,yAy^^ 



(C5) 



To go from the third to the fourth line here we use the 
completeness relation (|C3p . We have therefore shown 
that A 'P — M, which implies that the eigenvalues of M 
can only be zero or one. 

Our second step is to show that M is invertible, and 
that it can therefore have no zero eigenvalues. We do 
this by invoking the assumption of invertibility of the 
mapping from the {(?,,, p,;} to ipix), which implies that 
there exists some generalized matrix Cr/.x, of the same 
form as B,^ x, such that 



dx Cji 



St/j^x) 
Sip*{x) 



(C6) 



We then invoke our basic expression for Sipix) as a func- 
tion of the {(?,,, p^}, Eqn. (fT5)) . which may be compactly 
expressed using our generalized matrix A: 



dip{x) 
6^*{x) 



— Ax^ri 



(C7) 
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Combining these two expressions yields 



(C8) 



where we have again used completeness (|C3[) to go from 
the first line to the second, and in the last line we simply 
define 



= I dxC\ 



(C9) 



This final result holds for all possible qr],Pri, which implies 
that LM is the identity matrix. Hence the inverse of M 
exists (because we have constructed it, from the assumed 
C, as L). 

Since all eigenvalues of M are therefore either zero or 
one, and since none of them can in fact be zero, it fol- 
lows that they are all one. Thus M is the identity matrix, 
which is the orthonormality condition. This further im- 
plies that C = B. 



Appendix D: Derivation of the zero frequency 
solutions 



We present here in very detail how to obtain all 
(asymptotic) zero frequency solutions of the BdGE (|47)) 
and This is best achieved by rewriting the BdGE 

analogous to the factorization method of subsection IV CI 
Without loss of generality we set the soliton position xq 
to zero in this appendix. 

At first we solve the homogeneous equation (|Tf)) with 
the ansatz R{x) — e'^™^{f{x) + ig{x)) with real functions 
/ and g. With the definition of the operators Q and 
from (|43p . the second order differential equation in 
complex function space can be rewritten as a coupled 
system in real function space: 



"o' 










V2/3g 2(/32 



- K^) + 





7" 




9_ 



(Dl) 



Except for the imaginary unit the matrix involved is iden- 
tical to the one in pi)) . From the bottom line one finds 
that 



Qf = -V2Pg + ^ cosh^ 
v2 



(D2) 



where the second term on the right hand side is annihi- 
lated by for arbitrary constant C4. Inserting this in 
the top line of (|Dip it can be cast into an inhomogeneous 
second order differential equation for g: 



——dl — K^ scch^KX 



g — — K/3C4 cosh Kx 



(D3) 



If we define 

gi = tanh kx 
g4 = — cosh KX 

K 
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: tanh kx — 1 



(D4) 



one finds that gi, gs are the complementary functions 
of the differential equation (jD3l) and 54 is the particular 
solution. Hence, the general solution is given by 



(D5) 



where Ci, and C3 are arbitrary constants and C4 is the 
same solution parameter as in (|D2I) . 

With the general solution (IDSP one can search for a 
solution / in (|D2I) . The general solution to this inhomo- 
geneous first order equation is straightforwardly found to 
be 



4 

E 



(D6) 



where we have used the definitions 



/i 

U 



K 

P_ (I 
K \2 

4k2 



/2 = sech KX 



(tanh KX + KX sech^Kx) 



K 



(sinh 2kx + 3(tanh kx + kx sech^ra)) 

(D7) 



and introduced C2 as solution coefficient for the comple- 
mentary function /2 . This yields the four complementary 
solutions Rj — e~™^{fj + igj) (with 52 = 0) for the ho- 
mogeneous BdGE (|T7l) as presented in subsection I VI Al 
In the same manner one can transform the inhomo- 
geneous BdGE (gS]) with the ansatz S = e~™^(ig - p), 
where q and p are real functions, into 







. 9 _ 
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P. 



(D8) 



with the inhomogeneities given by (jD5|) and (|D6p . The 
lower line can be viewed as an inhomogeneous first order 
differential equation for the function Qq + \/2(3p. Its 
particular solution can be shown to be 



^Qq^-2f5p=^ + ^{ 



— sinh 2kx) 
2 ' 



Acos^\x 

K 



(D9) 



with the same constants C, as in (ID5I) and (ID6I). The 



homogeneous part of this equation is solved by g^, which 
would finally lead to a solution (0,ii?4) of the complete 
BdGE (gT]) and (gHl)- As we already know from IVI Al 
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that the sohitions of the type {0,iRj) exist, we do not 
consider it in the following, but concentrate on finding 
the particular solutions. 

Inserting (|D9[) into the upper line of (jPSp yields a sec- 
ond order differential equation for p of the form (|D3[) . 
but with different inhomogeneities: 



— — 9| — sech^Kxj p 
C2/2 - C3 (fs + ^ (kx + i sinh2Ka; 



C4 ( A —xcosh^Kx 



(DIO) 



Thus the complementary functions for the homogeneous 
part of this equation are exactly the same functions gi 
and gs, and the remaining problem is to determine the 
particular solutions for p. We decompose p = J2j=i ^jPj^ 
where the coefficients C, are chosen to be the same as in 



(|D5I) and (ID6p . and solve the equation (IDlOp for each 
individual j. Since (jPlOp does not contain an inhomo- 
geneity with Ci, the corresponding solution pi is zero, 
and for j — 2 one easily finds the particular solutions p2 , 
such that we have: 



Pi=0 



1 



P2 



(Dll) 



We have not found the exact solutions for p3 and p^. 
Nevertheless, one can expand (IDIOI) in powers of e^'^'^' to 
determine their asymptotic behavior for x — > ±00. This 
will be sufficient enough for our purpose since one can 
prove from their asymptotic form that they do not lead 
to physical zero modes, i.e. these modes are either not 
properly normalizable in the case of the infinite system 
without boundary conditions or they are incompatible 
with periodicity in case of the finite system. The leading 
terms are given by 



P3 
P4 



/3 



■ sgn xe 
^ + 3/32 



2k,\x\ 



16k4 



■ sgn x 



4^3 



(D12) 
(D13) 



and are due to the exponentially growing inhomogeneities 
in (ImOll . 

With the given (asymptotic) solutions pj one can solve 
the first order equation (|D9p in a straightforward manner 
with the ansatz q — ^jlj- The particular solutions 

for J = 1, 2 are exactly given by 



1 



'^^ = 2.2 



(tanh KX + KX sech^Kx) 
and for j = 3, 4 we find for x — > ±00 

'74 



2(3 
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(D14) 



16k4 



2N 

K X Sgn X 



(D16) 



to leading order in e^'''^'. Combining the results to Sj = 
Q-™x{iq. — Pj) yields the exact zero frequency solutions 
to (|48| with inhomogeneities Rj for j = 1,2, which are 
equivalent to the solutions presented in ([55)1 and ([55]) . In 
the case j — we correspondingly find the asymptotic 
solutions. It is important to recognize that we have found 
the asymptotic behavior of global solutions S3 and S'4. 



Appendix E: Linearization for a bright soliton 
background 

For attractive particle interaction the sign of the j^"!* 
term in ^ changes, and a time independent solution of 
the corresponding GPE in a frame moving with speed /3 



idt4>{x,t) 



]^^l-\^\''-^i-iP^,]^P{x,t) (El) 



is the bright soliton 

V;f,(a;) = e-*'^(^-^°^«;sechK(a; - a;o). (E2) 

with the chemical potential given by /i = — ('«^ + /3^)/2. 
The linearization around this bright soliton background 
analogous to Section HVl leads to the BdGE for the mode 
functions 



Hb = -Idl ~ 2\^b\'' - ^i - iPd, . 



(E3) 



Notice the sign changes in (|Eip and (|E3p due to the sign 
change of the particle interaction. Abbreviating x — x — 
xq the BdGE for unit mass have the non-zero frequency 
solutions 



Rri = {KtanhKxdx 



Sri — {KtSLnhKxdx H (fc ~ )) I 2 

2 k 



K^ sech^ Kx)d{x) 
d{i) 



(E4) 



d{x) 



sin kx 
cos kx 



where both sin kx and cos kx are possible for d{x) . The 
non-zero frequency spectrum is given by fife = k'^/2 + \fj,\, 
which reveals the macroscopic energy gap (k^ + (3"^)/ 2 
between the lowest lying zero energy soliton translational 
mode (of positive mass) 



(E5) 



(D15) 


Rbz 


_ 1 


'-{t/3 + dx)^:b 




Sbz 


~ V2k 


idp^b 



and the non-zero frequency modes (|E4I 
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